Coherent quantum ratchets driven by tunnel oscillations: Fluctuations and 

correlations 
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We study two capacitively coupled double quantum dots focusing on the regime in which one 
double dot is strongly biased, while no voltage is applied to the other. Then the latter experiences 
an effective driving force which induces a ratchet current, i.e., a dc current in the absence of a 
bias voltage. Its current noise is investigated with a quantum master equation in terms of the 
full-counting statistics. This reveals, that whenever the ratchet current is large, it also exhibits 
some features of a Poissonian process. By eliminating the drive circuit, we obtain a reduced master 
equation which provides analytical results for the Fano factor. 

PACS numbers: 73.23.Hk, 05.60.Gg, 72.70.+m 



I. INTRODUCTION 

The realization of double or triple quantum dots in 
a linear arrangement^^— or in a ring configurationA'^ 
enables transport experiments in which electrons flow 
through delocalized orbitals. This delocalization is 
visible in modified quadruple points of the charging 
diagramii"— Recently, also the capacitive coupling of two 
double quantum dots has been achieved.—"— If each dou- 
ble quantum dot is occupied by one electron, the setup 
represents two interacting charge qubits.— >^ Upon open- 
ing one double dot, a current flows and may be used for 
readout of the other double dot which still forms a charge 
qubit^-i^ 

Opening both double dots enables experiments with 
two interacting mesoscopic currents. For example, pre- 
dominantly coherently transported electrons perform 
tunnel oscillations which act on the other double dot as 
an effective ac force. This may induce a ratchet or pump 
effect, i.e., cause a net current in the absence of any bias 
voltage.ii A similar ratchet effect has been realized by 
coupling a double dot to a quantum point contactii This 
effect is closely related to Coulomb dra g^^d^ and to us- 
ing a double quantum dot as noise detectori^ d^d^ In turn, 
interacting channels may block each otherji^r— 

A common characterization of the current fluctua- 
tions in a mesoscopic conductor is the full-counting 
statistics^i^ whose cornerstone is a counting variable 
for the lead electrons. In this way, one obtains a 
cumulant-generating function for the transported charge. 
Of particular interest is its variance, because it relates to 
the zero-frequency limit of the current-current correla- 
tion functionJ^i Moreover, it indicates whether transport 
is sub or super Poissoniau)^ even though a more faithful 
criterion is the g^^-' functionj2^ Generally, some further 
cumulants are specific to the system while beyond a cer- 
tain order, cumulants exhibit universal featuresi^ 

In this work, we explore the noise properties of 
the ratchet mechanism proposed in Ref. 14 for capaci- 
tively coupled double quantum dots focusing on the full- 
counting statistics. Besides a numerical study with a 



master equation for the full ratchet-drive setup, we de- 
rive in the spirit of Ref. [1^ an effective master equation 
for the ratchet under the influence of the drive circuit. 
This provides an analytical expression for the cumulant 
generating function. This approach is beyond a more 
heuristic elimination of the drive circuit^^ and beyond 
a golden-rule calculation^ because it includes effects 
stemming from delocalization and from the broadening 
of the ratchet levels. Therefore, it holds also for small 
ratchet detuning. 

In Sec. ini we introduce our model Hamiltonian for 
the setup sketched in Fig. [T] and, moreover, introduce a 
quantum master equation for the full system. Section Hill 
is devoted to the elimination of the drive circuit which 
provides our analytical results. In Sec. IIVI we present 
our numerical results for the higher-order cumulants and 
test the quality of our approximations. 




FIG. 1. Quantum ratchet (lower circuit, unbiased) capac- 
itively coupled to a drive circuit (top) biased by a voltage 
Vb- Each circuit is modeled as two- level system with tunnel 
couplings Tra and Tdr, respectively. The ratchet possess a de- 
tuning e, while the drive circuit is undetuned. The dot-lead 
tunnel rates are Fra and Fdr, while /^a, a = 1, . . . , 4, denotes 
the chemical potentials of the leads. 
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II. MODEL AND METHOD 
A. Hamiltonian 

The four quantum dots and the leads (see Fig. [Ij are 
described by the Hamihonian H = Hs + Hb + Hv, where 



a<0 



(1) 



models the quantum dots with the electron creation and 
annihilation operators dj, and da- The ratchet circuit 
(a ~ 1,2) has inter-dot tunneling Tia and detuning e, 
such that the level splitting becomes 



5= Ve'+4|T,a|^ 



(2) 



The levels of the drive circuit (a = 3, 4) are not detuned 
and possess a tunnel matrix element T^r- The setup is 
assumed symmetric, such that inter-channel Coulomb re- 
pulsion reads U = U13 = C/24, while the internal re- 
pulsions U12 and U34 are assumed so large, that each 
channel can be occupied with at most one electron. The 
inter-channel coupling U by contrast, is relatively weak 
but nevertheless is the relevant interaction for inducing a 
ratchet current We do not take into account more indi- 
rect interaction mediated by phonons^^ or by a qubit.^*^ 
Each dot a is coupled to a lead with chemical potential 
/ia, where /ii — /i2, while /ia ^ /i4, such that all levels of 
the drive circuit lie within the voltage window. The lead 
Hamiltonian and the dot-lead couplings read 



ka 

Hv =51 + VkJlcko 



(3) 
(4) 



ka 



respectively, where cj^^and Cka a-i'e the fermionic opera- 
tors and Ska is the corresponding single-particle energy. 
The system-bath interaction Hy is determined by the ef- 
fective tunnel rates rQ(e) — ^nj^k \Vka\'^S{e—eka) = Ta, 
which within the wide-band limit are assumed energy- 
independent. Throughout this work we use units in which 



B. Cumulant generating function and master 
equation 

We are interested in the low- frequency properties of the 
transport process, which can be characterized by the dis- 
tribution of the number of transported electrons at large 
times or, equivalently, by the corresponding cumulants. 
We thus introduce for each lead a a counting variable 
Xa such that we obtain the moment generating function 
^gZx(JV-JVo)^^ = exp[G{x,t)], where x = (xi, X2, X3, X4), 



while N = {Ni, N2, N3, N4) is the electron number in 
each lead in vector notation with the initial value TVo = 
{N)t=o- Obviously, the Taylor coefficients of exp[G(x, t)] 
are the moments of the lead electron distributions. This 
allows the definition of occupation cumulants as Taylor 
coefficients of ln(e*'^'-'^~-'^''))f . Eventually they grow lin- 
early in time.^'^ Thus, the information about the station- 
ary limit is contained in the time derivative in the long- 
time limit, so that the (particle) current cumulants of 
leads 2 and 4 can be written as 



gn Q 



9(1x2)™ d{ixiY dt 



Gix,t) 



(5) 



Owing to charge conservation, the low-frequency prop- 
erties of the currents in leads 1 and 3 are identical with 
those of leads 2 and 4, respectively. Thus, it is suffi- 
cient to consider only the latter. The first-order current 
cumulants are the ratchet current = eoKi.o and the 
drive current I^,- = eoKo.i, where Bq denotes the elemen- 
tary charge. The second-order derivatives correspond to 
the zero-frequency limit of the current-current correla- 
tion functions,^* in particular, S^n = e^K2fi- Since our 
focus lies on the ratchet current, we introduce the nota- 
tion Kjn = Km,0- 

The cumulant generating function G{x,i) can be ob- 
tained from a Markovian master equation approach by 
unraveling the reduced density operator according to the 
number of electrons in the leadsi^ Alternatively, one 
may attribute a counting variable to the system-lead tun- 
nel operators^"— or multiply the full density operator by 
g^x -'v ijgfQpg tracing out the leadsi^ The resulting aug- 
mented density operator p(x,i) relates to the cumulant 
generating function via ln[tr /9(x. t)] = G'(x,i), while its 
limit p{x ^ is the usual reduced density operator. 
Moreover, p(x,t) obeys the master equation 



p[x.t) = C{x)p{x.t) 

a,s— ± 



^a PiX,t), (6) 



where Cq is the Liouville operator, in our case the one ob- 
tained within Bloch-Redfield approximation. For a short 
derivation, see Appendix |^ The superopcrator J'^ de- 
scribes tunneling of an electron from dot a to the respec- 
tive lead, while J'~ captures the opposite process. Thus, 
the counting variables keep track of the electron numbers 
in the leads despite that the latter are traced out. From 
the master equation ([5]), one can obtain numerically all 
cumulants within the recursive scheme of Ref. IsBI . Be- 
fore discussing these results, we aim at further analytical 
progress. 



III. ELIMINATION OF THE DRIVE CIRCUIT 

In order to reduce the number of degrees of free- 
dom, such that an analytically solvable master equation 
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emerges, we eliminate the drive circuit along the lines of 
Ref. [2^ We start by separating the master equation for 
piXit), Eq. ([6|), into contributions for the ratchet, the 
drive circuit, and their mutual interaction, 

P{X2,t) = ['Cra(X2)+^dr + C^^ra-dr]p(X2,i)- (7) 

Since we focus on the ratchet current, we keep only the 
counting variable X2 for the right lead of the ratchet cir- 
cuit. The interaction Liouvillian 



(8) 



is governed by the occupation imbalances Afiia = 712 — hi 
and Andr = hi — h^, which allow one to approxi- 
mately write the ratchet-drive interaction Hamiltonian 
asi^ U{hih^ -\- h2h/C) ~ {U /2) AndrAfira- Thereby we ne- 
glect terms that cause global shifts of all dot energies. 
They are not relevant here, because for all parameters 
considered below, the onsite energies stay far from the 
Fermi surfaces. 

After transforming master equation ([7]) into Laplace 
space, straightforward algebra^^ provides for the ratchet 
an effective Liouvillian Ccs, which follows from the rela- 
tion 

[z-C.^{X2,z)]-'^ir^.{[z-C{x2)]-'pT} (9) 

and depends on the stationary state p^^^^ of the drive 
circuit. Taylor expansion up to second order in the in- 
teraction constant U and subsequent evaluation of the 
partial trace yields 



/:eff(x2, z) = /:,a(x2) + UC^^^{z) + u^c^^iz) 

with the linear and quadratic corrections 



-(Andi.)[An„ 



X [An,a, .]l0[r^))((^^r^l[Anr 



(10) 

(11) 
(12) 



Here we employ the superoperator notation of Ref. l36l 
and define [M, •Jp = [Af, p]. Moreover, we have intro- 
duced the spectral decomposition of the ratchet Liouvil- 

with the eigenvalues A*-™' = 



lian, Em-^ra '\4>\ 
0, — Fra, — Fra/2 ± 1(5, and the left and right eigenvectors 



and l^^r"*)), respectively. Here a difficulty arises 
from the fact that the Liouvillian of a double quantum 
dot in the zero-bias limit is defective, i.e., it does not pos- 
sess a complete set of eigenvectors (for details, see Ap- 
pendix [B|). Then one may proceed either by constructing 
a generalized eigenbasis or by introducing a small pertur- 
bation that lifts the defectiveness, and finally consider the 
limit of vanishing perturbatiouj^ Since all levels are as- 
sumed to stay far from the Fermi surfaces, the impact of 
the interaction on the jump operators can be neglected 



safely. Thus, the jump operators J2 of the effective 
model coincide with the ones of the full Liouvillian. 

The dependence on the drive circuit enters via the 
Laplace transformed auto correlation function of the pop- 
ulation imbalance 



C{t) = (Andr(i)A«dr(0)) - (Andr) 



(13) 



evaluated at the eigenvalues of the ratchet Liouvillian, 
C{z - A^r^), which fulfills C*{z) = C{z*). For a deriva- 
tion, see Appendix [C] Below we will find that the poles 
of C{z) are related to the extrcma of the ratchet current. 

The linear contribution to the effective Liouvillian, 
L^^^ merely provides a small additional bias for the 
ratchet circuit, but does not induce any non-equilibrium 
effect. Thus, we omit this term and focus on the 

(2) 

impact of C)^^- By a straightforward but tedious 
evaluation of Eq. we obtain in the Fock basis 

1)(1|,|2)(2|,|2)(1|,|1)(2|} the expression 



/O 





A{z) 

VO B{z) 







B*{z*) 
A*{z*)J 



(14) 



where the terms 

(52 



A{z) 



s=± 



[2C{z + 

((5 4- sef 
4(52 



rra)-rraC'(z + F,.a)] 

C{z + T,.j2 + si5), 



(15) 



B{z) 



'(52 



[2C(z + r,.a)-rraC"(z + rra) 



- C{z + r,a/2 + ^(5) - C{z + r,j2 - iS)] (16) 



contain a non-Markovian correction through the depen- 
dence on the Laplace variable z. For the resulting effec- 
tive ratchet Liouvilhan, £ra(x2) + C^^'^lff (2;), the cumu- 
lant generating function can be obtained by a standard 
procedure, namely by computing the eigenvalue that van- 
ishes for X2 0^ This yields a somewhat bulky expres- 
sion and, thus, we restrict ourselves to the Markovian 
limit obtained by z — > 0. By differentiation with respect 
to X2, we obtain the current and the zero- frequency noise 
as 



eo26e 



Im[(-r,a/2-l-i(5)C(r,a/2-l-z(5)], (17) 



5ra ~ Ini [(-Frae' + tSs^ + td^)C{r,j2 + id)] , 



(18) 



where b = 4|ri.apC/^/^(4(52 + r^J. These expression, in 
turn, allow us to simplify the cumulant generating func- 
tion G{x2,t). After evaluating the time derivative and 
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FIG. 2. (a) Real and (b) imaginary part of the auto cor- 
relation function (dashed lines) of the drive population im- 
balance in Laplace representation, C{zg), evaluated at the 
broadened resonance of the ratchet Liouvillian, = Fra/S -f 
i{e^ + 4\Tr^\'^y/'^ with F^a = 0.5 Td^ and T^a = 0.2Tdr, as 
function of the detuning e. The dotted lines correspond to 
the limit Fra — >■ 0. The dot-lead coupling is Fdr ~ Tra- 
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FIG. 3. (a) Ratchet current Ira, its zero- frequency noise Sra, 
and (b) the resulting Fano factor F = Sra/eo|/ra| as func- 
tion of the tunnel matrix element of the drive circuit. The 
results are computed with the full master equation. The 
other parameters are Fra = Fdr = 0.1 Tra, U = 0.5 Tra, 
and e = 5 Tra- The vertical dotted line marks the resonance 
condition 4|Tdr|^ = e'^ + 4|Trap. 



the limit t ^ oo contained in definition ^ of the current- 
cumulant, we obtain the generating function 



»(J"ra/eo)sin(x2) + (5'i.a/eg)[cos(x2 
. . fcr2 



|^^[C0S(X2) - 1] 



(19) 



which within the present approximation contains the 
full information about the low-frequency properties of 
the ratchet current. The presence of the counting vari- 
able X2 in the denominator, however, renders the ac- 
tual calculation of higher-order cumulants a formidable 
task. Only in the golden-rule limit, i.e., to lowest order 
in Fra, the denominator becomes independent of X2, so 

that Gi^)(x2) = «(/r^a°V5o)sin(x2) + (5L°Veg)[cos(x2)-l]. 
Consequently, we obtain the current cumulants 



gn 



9(«X2 



for odd n 

X2=o I Sra' /sq for even n 

(20) 

where the upper index (0) refers to the limit Fra — > 0. 
It turns out to be a good approximation, unless uni- 
versal cumulant oscillations set in, as we will discuss in 
Sec. UVBl 

Before testing the quality of this approximation and 
the parameter dependence of the results, we close this 
section by a remark on a formal aspect of the perturba- 
tion theory. Both the current ([T71) and the zero-frequency 
noise (jlSp are proportional to the auto correlation func- 
tion of the population imbalance, C{z), evaluated at 
the broadened level splitting of the ratchet, where the 
Laplace variable reads Zg = V^a_/2 -|- i{e'^ + 4|ri.ap)^/^. 



Thus we expect the ratchet current to exhibit resonance 
peaks. Taking into account the broadening distinguishes 
the present result from that of Ref. 14. There the ratchet 
current has been computed from the golden-rule rates 
for noise- induced transitions between ratchet eigenstates. 
While this treatment accounts properly for delocaliza- 
tion effects, it predicts too pronounced resonance peaks. 
Formally, the golden-rule solution is restored by the re- 
placement C(rra/2 + i5) C{iS) in Eq. dTT]). Figured 
visualizes that for ratchet detunings close to resonances, 
the difference between the two approximations may be 
significant. 



IV. CHARACTERIZATION OF THE RATCHET 
CURRENT 

Before starting with the analysis of the ratchet current 
fluctuations, let us compare the present case to that of 
a ratchet driven by an ac field. There, the current ex- 
hibits resonance peaks with large current and low zero- 
frequency noise i^^i^ For large driving amplitudes, the 
same behavior is visible at multi-photon resonances. Fig- 
ure [3] shows the corresponding result for the present driv- 
ing by tunnel oscillations. When the level splitting of the 
ratchet matches the tunnel frequency of the drive circuit, 
we indeed observe the qualitatively same behavior. Here 
however, we do not find higher-order resonances and, 
moreover, the Fano factor does not reach the extremely 
small values found in Ref. [H. The reason for this is that 
for realistic parameters, the driving via Coulomb inter- 
action with the upper circuit is much weaker than direct 
ac driving by, e.g., a high-frequency gate voltage. '^^ The 
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Full master eq. Effective model . . .without broadening 
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FIG. 4. Ratchet current (upper row), zero-frequency noise (middle), and Fano factor (lower row) as function of the ratchet 
detuning e for various tunnel couplings Tra. Results for the full master equation (solid lines) are compared to the analytical 
results H17p and (|18p (dashed). Dotted lines mark the golden-rule results which ignore the broadening Fra. The dot-lead 
tunneling rates are Fra = 0.5 Tdr and Fdr = 0.2 Tdr, while the inter-channel coupling reads U = 0.2 Tdr. 



kink in the Fano factor stems from a small step in the 
current and can be attributed to an energy difference of 
many-particle states that crosses the Fermi level of the 
ratchet This confirms our picture in which the tunnel 
oscillations of electrons in the drive circuit act like an ac 
driving with (angular) frequency = |2Tdr| determined 
by the tunnel splitting. 



A. Zero-frequency noise and Fano factor 

If the tunnel coupling of the ratchet is smaller than 
that of the drive circuit, Tia < Tdr, one can adjust 
the ratchet bias e such, that the resonance condition 

+ 4|Trap = 4|TdrP is met. By contrast, for T,^ > Tdr 
this is not the case. In order to first sketch the global be- 
havior, we first consider the current, the zero-frequency 
noise, and the resulting Fano factor in dependence of 
the ratchet bias. We compare numerical results obtained 
from the full master equation with the analytical solu- 
tion of Sec. mil Moreover, we also discuss the analytical 



expressions p7)) and ^TE\\ to lowest order in Fra, because 
this restores the golden-rule results of Ref. 

Figure |4] provides an overview to the behavior. The 
current which is depicted in the first row, exhibits the 
expected resonance peaks provided that Tra < Tdr- If 
the tunnel matrix element is rather small (Tra = 0.2 Tdr), 
we witness also the small peaks at small values of e, which 
we predicted within our analytical treatment. Upon in- 
creasing the inter-dot tunneling Tra, the current peaks 
naturally increase as well. Once Tra > Tdr, the reso- 
nance peaks fade away while the structure at e ~ be- 
comes rather pronounced. In all regimes, the analytical 
result (|17p for the current is well confirmed. The main 
difference is the absence of the slight asymmetry with re- 
spect to reverting the detuning, e — )• —e. Nevertheless, 
the characteristics of the current as function of e is by 
and large antisymmetric, which implies a current rever- 
sal close to zero detuning. In order to capture also the 
lack of perfect antisymmetry, we would have to consider 
the linear perturbation C^^fj which, however, would im- 
pede concise analytical results. For the parameters used 
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FIG. 5. Cumulant ratio versus order n for the pa- 

rameters used in the second column of Fig. |4] for two values 
of the detuning e. The value e ~ 1.94 Tdr corresponds to the 
resonance between ratchet and drive circuit. The horizontal 
lines in panel (a) mark the analytical result (|20p valid to low- 
est order in Fra, i.e., F and 1/F. The dotted lines serve as 
guide to the eye. 



in Fig. the golden-rule expression of Ref. [T3, i.e., 
Eq. P?|) to lowest order in Fia, reproduces the behavior 
only qualitatively. It predicts too sharp peaks, because 
this approximation does not account for the level broad- 
ening of the ratchet. The deviation is quite significant in 
the non-resonant case T^a > Idr. 

The main features such as the location of the peaks 
are also found for the zero-frequency noise plotted in 
the middle row. An important difference is found only 
close to £ = 0, where the current vanishes for symmetry 
reasons. The noise nevertheless remains finite and may 
even have a peak. This behavior is reflected by the Fano 
factor which stays close to the Poissonian value F — 1 
for detunings far from the current reversal point £ = 0. 
There the current vanishes, while the noise remains fi- 
nite, such that F diverges. The reason for this universal 
behavior can be understood from the analytical results 
for the current and the noise. Both and S'j-a depend on 
the drive circuit via the drive correlation function C{z) 
which, thus, cancels in their ratio, i.e., in the Fano factor. 
On a smaller scale, we observe in the Fano factor occa- 
sional kinks in less important regions in which the current 
is rather small. There a small change in the denominator 
of F = S'ra/eol^ral may have a strong effect. 
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FIG. 6. Correlation coefficient r versus detuning e for the 
parameters used in Fig. O 

B. Higher-order cumulants 

For a refined study of the current noise, we investigate 
also the cumulants of higher order, where we express the 
results in terms of the ratio between subsequent cumu- 
lants, For this quantity, the limit of small Fra 
is rather interesting, because our analytical result ([20)1 
implies that the cumulant ratio alternates between the 
Fano factor and its reciprocal. Such behavior is charac- 
teristic for a bi-directional Poisson process.'*'^ The usual 
Poisson process with |k„+i/k„| = const emerges as spe- 
cial case if the backward current is negligible. The higher- 
order cumulants of the effective Liouvillian for larger val- 
ues of Fra can be evaluated from the generating function 
(Uni), but the expressions become rather bulky, so that 
one has to resort to a numerical evaluation. 

Figure [5] shows a comparison of these two approxima- 
tions together with the result of the full master equa- 
tion. For a small ratchet detuning below the resonance 
[panel(a)], we find that all three solutions agree quite 
well and that the first few cumulants exhibit the pre- 
dicted alternation between the values F and 1/F. For 
higher orders, the generic universal cumulant oscillations 
set which obviously is beyond our analytical ap- 

proach. At the resonance, the universal oscillations start 
even already at lower order and Eq. (l20l) no longer holds. 
This is in agreement with our earlier observation that the 
broadening of the ratchet levels plays a significant role for 
resonant driving. 



C. Cross correlations 

Let us finally consider the cross correlations between 
the drive current and the ratchet current. They can be 
characterized by the cumulant ku, which is equivalent 
to the covariance of the transported charge in the two cir- 
cuits. As a dimensionless measure, we introduce the cor- 
relation coefficient r = kii/^K2o '^02, which is bounded 
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by — 1 < r < 1. The results depicted in Fig. |6] demon- 
strate that the correlation between the two currents is 
rather low. While it can be up to \r\ ^ 0.1 at the res- 
onances, it is hardly noticeable in the non-resonant case 



V. CONCLUSIONS 

A double quantum dot with detuned energy levels may 
act as quantum ratchet or quantum pump when driven 
out of equilibrium. An external force with zero net bias 
acting locally upon such system can induce a dc cur- 
rent. Here we investigated a ratchet with a particular 
driving, namely one that stems from the capacitive cou- 
pling to a further double quantum dot which, however, is 
biased. Electrons flowing through the drive circuit per- 
form tunnel oscillations which indeed induce phenomena 
similar to those induced by deterministic ac driving. In 
this work we mainly focused on the fluctuations of the 
emerging ratchet current. Besides a numerical solution 
with a master equation for all four quantum dots, we de- 
rived an effective ratchet Liouvillian by eliminating the 
drive circuit. In this way, we obtained analytical results 
even for higher-order cumulants, which agree well with 
those of the full master equation provided that the tun- 
nel splitting of the drive circuit is larger than that of the 
ratchet. 

As a common feature of driving by tunnel oscillations 
and driving by an ac field, we found resonance peaks at 
which the current assumes a maximum, while the rela- 
tive noise characterized by the Fano factor is minimal. 
However, clearly sub-Poissonian noise is only found for 
large detuning of the ratchet levels. This noise reduc- 
tion should be measurable, even though it is not as pro- 
nounced as in the case of ac driving, mainly because it re- 
quires large driving amplitudes which cannot be achieved 
by capacitive coupling. For less detuned ratchet levels, 
the Fano factor is typically of the order one, unless the 
detuning is so small that its orbitals are fully delocalized. 
Then the lack of sufficiently strong asymmetry keeps the 
current at a low value, while the zero-frequency noise 
stays finite. Thus, the Fano factor being the ratio of these 
two quantities assumes very large values. This generic 
behavior of the Fano factor is explained by our analyt- 
ical results which reveal that both the current and the 
zero-frequency noise are proportional to the correlation 
function of the drive circuit. Thus the Fano factor de- 
pends only on the shape of the ratchet eigenfunctions, 
while the correlation function cancels. 

The higher-order cumulants tend to alternate between 
two values. This indicates a bi-directional Poisson pro- 
cess and implies that a backward current becomes rel- 
evant. With increasing order, however, universal oscil- 
lations with ever larger amplitude dominate. The onset 
of the universal oscillations marks the point at which 
our analytically obtained higher-order cumulants signif- 
icantly deviate from those for the full master equation. 



Nevertheless, the physically relevant cumulants of lower 
order are well within our analytical treatment. 

The more global picture is such that the noise is close 
to the Poissonian level, whenever the current is relatively 
large. Thus possible applications and measurements of a 
ratchet current induced by tunnel oscillations, should not 
be hindered by current fluctuations. 
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Appendix A: Liouvillian and jump operators 

For a system-bath Hamiltonian, one can derive for the 
reduced system density operator p the Bloch-Redfield or 
Born-Markov master equation^"— 

p{t) ^ - i[Hs, p{t)] 

- / dsiiB[Hv, [Hvhs), p{t)Ro]] (Al) 
Jo 

= Cop{t), 

where Rq is the reference density operator of the en- 
vironment, while Hv describes the system-bath cou- 
pling. It has to be weak, such that coherent system 
dynamics dominates. Augmenting the environment den- 
sity operator by a counting variable for the lead elec- 
trons according to p{t)RQ p{t)e'^^^'^~-'^°^ Rq, yields 
the x-dependent density operator whose trace is 

the moment-generating function introduced in Sec. IIIII — 
Furthermore, one obtains by inserting the same ansatz 
into Eq. (|A1[) the master equation ^ upon which all our 
results are based. 

For the evaluation of the time integral in Eq. (jAll) , it is 
convenient to work in the eigenbasis of the system Hamil- 
tonian, defined by Ek\m) — Hs\m). Then one obtains 
for the density matrix elements the equation of motion 
Pmn(x,0 = J^kii'^Mlmn.ki Pkiix,^), whcrc the decom- 
posed Liouvillian 

[^ix)]mn,kl — [^o]mn,kl + ^ (e**^° — l)[J^]mn,fci 

a.s— it 

(A2) 
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consists of the contributions 



Q.S — ± 

pk 



pn 



(A3) 



[Jq ]mn,kl — ■^i'Ja,mk + "fa,nl)[dl,]mk[da]ln , (A4) 
[Ja]m7i,kl — 2(70, fern + 7a,/n) [^a] mfc MqI/k ■ (A5) 



The effective rates 



(A6) 
(A7) 



depend on the Fermi functions of the leads, fa{ijj) 
f{uj - fia) = {cxp[(a; - fia)/kBT] + while 



^2n J2\Vkc.\^S{u:-Skc.) 



(A8) 



denotes the spectral densities of the dot-lead couplings, 
which within a wide-band limit are assumed energy in- 
dependent. 



Appendix B: Spectral decomposition of the ratchet 
Liouvillian 

In the energy basis {|0)(0|, |e)(e|, \g){gl \g){e\, \e){g\} the 
ratchet Liouvillian for vanishing counting variable X2 
reads 



r — 





r 

^ ra 











\ 





_r 

J- ra 






































2 ^ 


-iS 




V 











2 





(Bl) 



Within the perturbative treatment of Sec. IIIIl we need 
to compute functions /(£ia) of this matrix, such as the 
propagator exp(£rai) or the resolvent (z — £ra)~^, which 
is usually achieved by spectral decomposition of the Li- 
ouvillian. Here however this is hindered by the fact that 
£ra is defective, i.e., it does not possess a complete set of 
eigenvectors. The problem arises from the upper block 



(B2) 





Tra 





^ ra 


Tra 






which we transform via 



5': 



^0 -1 

~i/r,, 
,1 1 i/r,.a 



(B3) 



to the Jordan canonical form^i^ 



(B4) 



Its eigenvalues obviously are and the twofold degenerate 
— Fra, and one immediately finds two vectors that obey 
the eigenvalue equation, namely 



i|0>=0|0), 
ill) = -F,, 



ID 



(B5) 
(B6) 



A generalized eigenbasis can be found by including a 
third vector |2) that fumils^"^ 



i|2) = -Fra|2) + |l), 



(B7) 



i.e., one adds the eigenvector of the degenerate sub- 
space. By repeated multiplication with L follows i'"' 1 2) = 
(-F,a)'=|2) +/c(-F,a)'^'-i|l), which implies 



/(L)|2)=/(A2)|2)+/'(A2)|1) 



(B8) 



where both the function of a matrix and its derivative are 
defined as the corresponding Taylor series. This relation 
together with the usual f{L)\k) — /(Afe)|fc) for fc = 0, 1, 
allows us to evaluate any ,f{L). In particular, we find the 
propagator 



and the resolvent 



e-r-* 
-i^-* 1 - (1 + F,.at)e-r-* 



(B9) 







(BIO) 



A poor man's approach to this procedure^l is to in- 
troduce a small perturbation that lifts the degeneracy 
of L. After evaluating /(L), one considers the limit of 
vanishing perturbation. 



Appendix C: Stationary state of the drive circuit 

The effective ratchet Liouvillian derived in Sec. IIIIl 
results from a perturbation theory with the stationary 
states representing the zeroth order. They are deter- 
mined by the master equation^ 

Pdr = -i[^S,dr, Pdr] + 'D{d\)pdt + T>{d4)pdr, (CI) 



9 



with the system Hamiltonian of the drive 



ance Anji- = 0,4 — 713 in Laplace space is defined as 



(C2) 



The last two terms in the master equation describe dot- 
lead tunneling, which in the limit of a large voltage bias 
obeys the Lindblad form 

'D{d4)pdr = Tdr d4Pdr4 ~ -^{4«^4, Pdr} ■ (C3) 

For the drive circuit, the stationary state can be obtained 
conveniently after a decomposition of the Liouvillian into 
the corresponding Fock basis, by which we find 

Stat _ }_ 

4|TdrP 

XI r2^ + 4|Tdr|2 -2irdrTd*, I . (C4) 
2iTdrTd, 4|Td,|2 

The auto correlation function of the population imbal- 



C{z)^ / dte-^*[(Andr(t)Andr(0))-(Andr)'] 
Jo 

= (Andr(2)Andr) - -(Andr)^ (C5) 



with the stationary occupation 



(Andr) 



^ dr 



+ 12|TdrP 

and the corresponding correlation function 

1 2z + Fdr 



zTl + U\Tdr\ 



(C6) 



^ rg,(^ + rdr)^ + 4z(2^ + 3rd.)|rdrP . . 

(z + rd,.)2(2z + Fdr) + 4(2z + 3rdr)|TdrP ' ^ ' 
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